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[Name of Document] CLAIMS 
[Claim 1] 

An electromagnetic field analyzer comprising: 

dividing means for dividing form data as an analysis 
object into coarse elements and fine elements; 

forming means for forming a matrix making an 
electromagnetic field vector of the coarse elements divided 
by the dividing means related to an electromagnetic field 
vector of the fine elements, and 

calculating means for calculating an approximate 
solution of the electromagnetic field vector of the fine 
elements by applying an iteration method of simultaneous 
linear equations while referring to the matrix formed by the 
forming means. 
[Claim 2] 

An electromagnetic field analyzer as claimed in claim 1, 
wherein the forming means forms the matrix by expressing the 
elements of electromagnetic field vectors at the sides at a 
fine element using an interpolation function in the coarse 
elements . 
[Claim 3] 

An electromagnetic field analyzer as claimed in claim 2, 
wherein, when the length of a side 1^ of the fine element is 
|lj(. an interpolation function showing the relation between 
the electromagnetic field at position x of the side 1^ of the 
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fine element and the electromagnetic field at a side j of 
the coarse element is n/ (x) , and the tangential vector of 
the side 1^ of the fine element is t^, the forming means 
forms a matrix P^^ by using the following expression. 
Expression 1 



[Claim 4] 

An electromagnetic field analyzer as claimed in any one 
of claims 1 to 3, wherein, in the calculation means, the 
accuracy of an approximate solution of the electromagnetic 
field vector of the fine elements is improved such that 
high-frequency components included in the approximate 
solution of the electromagnetic field vector of the fine 
element are removed by applying a stationary iteration 
method of simultaneous linear equations, that a residual in 
the fine elements is mapped to a residual in the coarse 
elements by using the matrix formed by the forming means, 
that a correction vector to the coarse elements is formed by 
applying a direct method or a non • stationary iteration 
method of simultaneous linear equations, and that a 
correction vector to the fine elements is obtained by using 
the matrix formed by the forming means. 
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(Claim 5] 

A computer program to cause a computer to execute an 
electromagnetic field analyzing method for analyzing an 
electromagnetic field as an analysis object, 

wherein the computer contains a first storage means for 
storing divided elements and a second storage means for 
storing a matrix, and 

wherein the electromagnetic field analyzing method 
contains the steps of dividing form data as an analysis 
object into coarse elements and fine elements and storing 
the elements in the first storage means, forming a matrix 
making an electromagnetic field vector of the coarse 
elements stored in the first storage means related to an 
electromagnetic field vector of the fine elements and 
storing the matrix in the second storage means, and 
calculating an approximate solution of the electromagnetic 
field vector of the fine elements by applying an iteration 
method of simultaneous linear equations while referring to 
the matrix stored in the second storage means. 
[Claim 6] 

A storage medium, being computer readable, for 
recording a program to cause a computer to execute an 
electromagnetic field analyzing method for analyzing an 
electromagnetic field as an analyzing object, 

wherein the computer contains a first storage means for 
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storing divided elements and a second storage means for 
storing a matrix, and 

wherein the electromagnetic £ield analyzing method 
contains the steps of dividing form data as an analysis 
object into coarse elements and fine elements and storing 
the elements in the first storage means, forming a matrix 
making an electromagnetic field vector of the coarse 
elements stored in the first storage means related to an 
electromagnetic field vector of the fine elements and 
storing the matrix in the second storage means, and 
calculating an approximate solution of the electromagnetic 
field vector of the fine elements by applying an iteration 
method of simultaneous linear equations while referring to 
the matrix stored in the second storage means. 



f 
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(Name of Document] SPECIFICATION 

[Title of the Invention] ELECTROMAGNETIC FIELD ANALYZER* 
ELECTROMAGNETIC FIELD ANALYZING PROGRAM, AND STORAGE 
MEDIUM FOR RECORDING THE PROGRAM 

[Technical Field] 
[0001] 

The present invention relates to electromagnetic field 
simulation technology used in development and designing of 
high • frequency products such as multilayer products, giga 
filters, EMI (Electro*Magnetic Interference) elimination 
filters, etc., and more particularly to an electromagnetic 
field analyzer, an electromagnetic field analyzing program, 
and a storage medium for recording the program in which the 
simulation time can be shortened. 
[Bac]cground Art] 

[0002] 

In recent years, the development and designing of high- 
frequency products have been vigorously performed and, in 
order to perform an electromagnetic field analysis of high- 
frequency products, electromagnetic field simulation 
softwares have been widely used. Generally, there are many 
cases where a finite element method is used for 
electromagnetic field analysis. In the electromagnetic 
field simulation using the finite element method, since it 
takes a long time to solve simultaneous linear equations. 
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high-speed solving methods are desired. In particular, in 
the field of electromagnetic field analysis, since nothing 
except for a direct method called Gauss' elimination method 
can be used, the problem of taking a long calculation time 
is more serious. 



However, by a paper (R. Hiptmair, "Multigrid method for 
Maxwell's equations." SIAM Journal of Numerical Analysis, 
vol. 36, no. 1, pp. 204-225, 1999) made public by Hiptmair 
in 1999, it is proved that a solution method called a 
multigrid method can be also used in an electromagnetic 
field analysis and the method is paid attention to as a 
greatly speeding up solution when compared with a direct 
method . 



In the multigrid method, an analysis object is divided 
into two classes of fine elements and coarse elements and 
the solution of fine elements is found by making use of the 
solution of coarse elements. 



When simultaneous linear equations are calculated by 
using Gauss' elimination method, generally the calculation 
time is proportional to third power of dimensions of a 
matrix. Since the size of dimensions corresponds to the 
number of elements, when the number of elements becomes 



[0003] 



[0004] 



[0005] 
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double, the calculation time becomes 8 times as large, and, 
when the number of elements becomes 10 times as large, the 
calculation time becomes 1000 times as large. Then, when 
Gauss' direct method is applied to only coarse elements, an 
approximate solution to fine elements is obtained by using a 
solution to coarse elements, and Gauss' elimination method 
is not used, the calculation time for fine elements can be 
reduced to 1/8 and 1/1000, respectively. Since additional 
calculation time is added, although the practical 
calculation time is not reduced as such, the calculation 
time can be greatly reduced without fail. 
(00061 

The multigrid method is classified into a method using 
a nested mesh and a method using a non-nested mesh. However, 
since a prolongation matrix is incomplete as described later, 
there is nothing except for using a method using nested mesh. 

[0007] 

Fig. 12 is drawings for describing the division into 
elements using a nested mesh. Fig. 12 (a) shows the 
division into elements using a nested mesh in the case of 
two-dimensional analysis where a coarse triangular element 
is uniformly divided to form fine triangular elements. 
Furthermore, Fig. 12 (b) shows the division into elements 
using a nested mesh in the case of three-dimensional 
analysis where a coarse tetrahedral element is uniformly 
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divided to form eight fine tetrahedral elements. As is 
understood from Fig. 12, there are geometrical restrictions 
between coarse elements and fine elements in the case of 
nested mesh. 

[Non-Patent Document II R. Hiptmair, "Multigrid method 
for Maxwell's equations," SIAM Journal of Numerical Analysis, 
vol. 36, No. 1, pp. 204-225, 1999. 

[Non-Patent Document 2] D. Dibben and T. Yamada, "Non- 
nested multigrid and automatic mesh coarsening for high 
frequency electromagnetic problems** lEEJ Investigating 
Research Committee material, SA-02-34, pp. 71-75, 2002. 
(Disclosure of the Invention] 

[Problems to be Solved by the Invention] 

[0008] 

Fig. 13 is drawings for describing a first problem of 
nested mesh. Fig, 13 (a) shows how a circle is divided into 
elements using a nested mesh to form coarse elements. 
Furthermore, Fig. 13 (b) shows how coarse elements are 
divided into fine elements using a nested mesh. As is 
understood from (a) and (b) in Fig. 13, a curved surface 
cannot be correctly expressed because of geometrical 
restrictions of a nested mesh. 

[0009] 

Furthermore, Fig. 13 (c) shows how a circle is divided 
into elements using a non-nested mesh to form fine elements. 





- 5 • 



JP Application No. 2003-412074 

As is understood from Fig. 13 (c) , since there are no 
geometrical restric tions « a curved surface can be correctly 
expressed . 



Fig. 14 is to describe a second problem of nested mesh. 
Fig. 14 (a) shows how a square is divided into elements by 
using a nested mesh to form coarse elements. Furthermore, 
Fig. 14 (b) shows how coarse elements are divided into fine 
elements. As is understood from (a) and (b) in Fig. 14, 
coarse elements cannot be partially divided into fine 
elements because of the geometrical restrictions of a nested 
mesh . 



Furthermore, Fig. 14 (c) shows how a square is divided 
into elements by using a non-nested mesh to form fine 
elements. As is understood from Fig. 14 (c), coarse 
elements can be partially divided into fine elements because 
of no geometrical restrictions. 



Host industrial products have curved surfaces as in a 
circular cylinder and a square and they inevitably contain 
partially fine portions. Accordingly, it is difficult to 
apply a nested mesh to industrial products. On the other 
hand, since there is no geometrical restrictions in a non- 
nested mesh, fine elements can be freely formed. However, 



(0010] 



[0011] 



[0012] 
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since the positional rellition is not systematic, it is 
difficult to make the electromagnetic field of coarse 
elements related to the electromagnetic field of fine 
elements. It is a prolongation matrix to make coarse 
elements related to fine elements, but no precise 
prolongation matrix has been found. This means that a 
multigrid method cannot be practically used, and accordingly, 
nothing can be used except for a Gauss' direct method which 
takes much time in the analysis of magnetic field. 
[0013] 

The present invention has been done to solve the above 
problem, and it is an object of the present invention to 
provide an electromagnetic field analyzer, an 
electromagnetic field analyzing program, and a recording 
medium having the program recorded therein, in which the 
analysis of an electromagnetic field can be made by using a 
multigrid method using a non-nested mesh. 

[0014] 

It is another object of the present invention to 
provide an electromagnetic field analyzer, an 
electromagnetic field analyzing program, and a recording 
medium in which the program has been recorded, in which the 
time required for an electromagnetic field analysis is 
greatly reduced. 

(Means for Solving the Problems] 



m 
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[0015] 

According to an aspect of the present invention, an 
electromagnetic field analyzer comprises dividing means for 
dividing form data as an analysis object into coarse 
elements and fine elements; forming means for forming a 
matrix making an electromagnetic field vector of the coarse 
elements divided by the dividing means related to an 
electromagnetic field vector of the fine elements, and 
calculating means for calculating an approximate solution of 
the electromagnetic field vector of the .fine elements by 
applying an iteration method of simultaneous linear 
equations while referring to the matrix formed by the 
forming means . 

[0016] 

Preferably, the forming means forms the matrix by 
expressing the elements of electromagnetic field vectors at 
the sides at a fine element using an interpolation function 
in the coarse elements. 

[0017] 

More preferably, when the length of a side 1^ of the 
fine element is |1J, an interpolation function showing the 
relation between the electromagnetic field at position xof 
the side 1^ of the fine element and the electromagnetic field 
at a side j of the coarse element is (x) , and the 
tangential vector of the side 1^ of the fine element is t|, 
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the forming means (25) £orms a matrix P^^ by Expression (14) 
to be described later. 
[0018] 

More preferably, in the calculation means, the accuracy 
of an approximate solution of the electromagnetic field 
vector of the fine elements is improved such that high- 
frequency components included in the approximate solution of 
the electromagnetic field vector of the fine element are 
removed by applying a stationary iteration method of 
simultaneous linear equations, that a residual in the fine 
elements is mapped to a residual in the coarse elements by 
using the matrix formed by the forming means, that a 
correction vector to the coarse elements is formed by 
applying a direct method or a non- s tat ionary iteration 
method of simultaneous linear equations, and that a 
correction vector to the fine elements is obtained by using 
the matrix formed by the forming means. 

[0019] 

According to another aspect of the present invention, 
in a computer program to cause a computer to execute an 
electromagnetic field analyzing method for analyzing an 
electromagnetic field as an analysis object, the computer 
contains a first storage means for storing divided elements 
and a second storage means for storing a matrix, and the 
electromagnetic field analyzing method contains the steps of 
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dividing form data as an analysis object into coarse 
elements and fine elements and storing the elements in the 
first storage means, forming a matrix making an 
electromagnetic field vector of the coarse elements stored 
in the first storage means related to an electromagnetic 
field vector of the fine elements and storing the matrix in 
the second storage means, and calculating an approximate 
solution of the electromagnetic field vector of the fine 
elements by applying an iteration method of simultaneous 
linear equations while referring to the matrix stored in the 
second storage means. 
10020] 

According to another aspect of the present invention, 
in a storage medium, being computer readable, for recording 
a program to cause a computer to execute an electromagnetic 
field analyzing method for analyzing an electromagnetic 
field as an analyzing object, the computer contains a first 
storage means for storing divided elements and a second 
storage means for storing a matrix, and the electromagnetic 
field analyzing method contains the steps of dividing form 
data as an analysis object into coarse elements and fine 
elements and storing the elements in the first storage means, 
forming a matrix making an electromagnetic field vector of 
the coarse elements stored in the first storage means 
related to an electromagnetic field vector of the fine 
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elements and storing the matrix in the second storage means, 
and calculating an approximate solution of the 
electromagnetic field vector of the fine elements by 
applying an iteration method of simultaneous linear 
equations while referring to the matrix stored in the second 
s torage means . 
[Advantages ] 
[0021] 

According to an aspect of the present invention, since 
the forming means forms a matrix making an electromagnetic 
field vector of coarse elements divided by the dividing 
means related to an electromagnetic field vectar af fine 
elements, it became possible to perform an electromagnetic 
field analysis by using a multigrid method using a non- 
nested mesh. Furthermore, since the calculating means 
calculates an approximate solution of the electromagnetic 
field vector of fine elements while referring to the matrix 
formed by the forming means, it is not required to calculate 
the electromagnetic field of fine elements by using a direct 
method, it became possible to largely reduce a time needed 
for an electromagnetic field analysis. 

[0022] 

Furthermore, since the forming means forms a matrix 
such that an electromagnetic field vector element at a side 
of fine elements is expressed by using an interpolation 





- 11 • 



JP Application No. 2003-412074 

function in coarse elements, it became possible to easily 
form a matrix. 



Furthermore, in the forming means, since a matrix P^^ is 
formed by Expression (14) to be described later, a matrix 
can be more easily formed. 



Furthermore, in the calculating means, an approximate 
solution of an electromagnetic field vector of fine elements 
is corrected by using a matrix formed by the forming means, 
it became possible to improve the accuracy of the 
approximate solution of an electromagnetic field vector of 
fine elements at high speed. 
(Best Mode for Carrying Out the Invention] 



Fig. 1 is an external appearance of an electromagnetic 
field analyzer according to an embodiment of the present 
invention. The electromagnetic field analyzer contains a 
computer main body 1, a display device 2, an FD (Flexible 
Disk) drive 3 to which an FD 4 is attached, a keyboard 5, a 
mouse 6, a CD-ROM (Compact Disk- Read Only Memory) drive 7 
to which a CD-RON 8 is attached, and a network 
communications device 9. An electromagnetic field analyzing 
program is supplied through a recording medium such as an FD 
4, a CD-ROM 8, etc. An electromagnetic field analysis is 



[00231 



[0024] 



[0025] 
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performed such that the electromagnetic field analyzing 
program is executed by the computer main body 1. 
Furthermore, the electromagnetic field analyzing program may 
be supplied to the computer main body 1 by anther computer 
through the network communications device 9. 
[0026] 

Fig. 2 is a bock diagram showing the construction of 
the electromagnetic field analyzer of an embodiment of the 
present invention. The computer main body 1 shown in Fig. 1 
contains a CPU (Central Processing Unit) 10, a ROM (Read 
Only Memory) 11, a RAM (Random Access Memory) 12, and a hard 
disk 13. The CPU 10 conducts processing while the CPU 10 
inputs and outputs data to the display device 2, the FD 
drive 3, the keyboard 5, the mouse 6, the CD-ROM device 7, 
the network communications device 9, the ROM 11, the RAM 12, 
or the hard disk 13. The electromagnetic field analyzing 
program in the FD 4 or the CD-ROM 8 is stored in the hard 
disk 13 through the FD drive 3 or the CD-ROM 7. The 
electromagnetic field analysis is performed such that the 
CPU 10 properly loads the electromagnetic field analyzing 
program from the hard disk 13 in the RAM 12 and executes the 
program . 

[0027] 

Fig. 3 is a block diagram showing the functional 
structure of the electromagnetic field analyzer according to 




- 13 - 

JP Application No. 2003-412074 

an embodiment of the present invention. The electromagnetic 
field analyzer contains a form data storage portion 21 for 
storing form data as an analysis object, a division- into- 
element portion 22 for dividing the form data such that the 
size of elements is changed, an element storage portion 23 
for storing the elements divided by the division - into- 
element portion 22, a vec tor /ma t r ix forming portion 24 for 
forming a vector and matrix in simultaneous linear equations, 
a prolongation matrix forming portion 25 for forming a 
prolongation matrix, a vector/matrix storage portion 26 for 
storing the vector and matrix formed by the vector/matrix 
forming portion 24, a prolongation matrix storage portion 27 
for storing the prolongation matrix formed by the 
prolongation matrix forming portion 25, an approximate value 
calculation portion 28 for calculating an approximate value 
of an electromagnetic field vector, an approximate solution 
correction portion 29 for correcting the approximate 
solution calculated by the approximate value calculation 
portion 28, and an electromagnetic • field vector storage 
portion 30 for storing as an electromagnetic field vector 
the approximate solution after having been corrected by the 
approximate solution correction portion 29. 
[0028] 

Moreover, the form data storage portion 21, the element 
storage portion 23, the vector/matrix storage portion 26, 
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the prolongation matrix storage portion 27, and the 
electromagnetic £ield vector storage portion 30 are provided 
in a fixed area in the RAM 12 or the hard disk 13 shown in 
Pig. 2. Furthermore, the function of the division- into- 
element portion 22, the vec tor /ma trix forming portion 24, 
the prolongation matrix forming portion 25, the approximate 
value calculation portion 28, and the approximate solution 
correction portion 29 is realized such that the CPU 10 shown 
in Fig. 2 executes the electromagnetic field analyzing 
program loaded into the RAH 12. 
[0029] 

Fig. 4 is a flow chart for describing the procedure of 
the electromagnetic field analyzer according to an 
embodiment of the present invention. First, the division- 
into-element portion 22 reads form data as an analysis 
object stored in the form data storage portion 21, divides 
the form data into elements, and preserves the divided 
elements in the element storage portion 23 (Sll) . The 
division - into * element portion 22 divides the form data into 
triangular or quadrangular elements in the case of two- 
dimensional analysis, and divides the form data into 
tetrahedral , triangle -pole - shaped, or square -pole • shaped 
elements in the case of three-dimensional analysis. 

[0030] 

The division - in to - element portion 22 divides the form 
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data while changing the size of elements. In order to make 
the description simple, here elements of two different sizes 
are formed and it is assumed that a set of larger (coarse) 
elements is made and a set of smaller (fine) elements is 
made O^. Moreover, a superscript of c represents coarse, 
and a superscription of F represents fine. 
[00311 

The sets of C and may be independently formed, but, 
when the set of O"" is first formed and then the set of is 
formed by subdividing the elements, the processing time is 
shortened. Regarding the method of subdivision, there are a 
method for subdividing elements which are larger in volume 
and a method for analyzing elements by using 0° and 
subdividing elements having large error. Either of them may 
be used, but the method in which error is evaluated has a 
better calculation accuracy. Moreover, regarding the method 
for evaluating error. Document No. 1 (Norio Kamiya et al . , 
"computer and application's monk, special edition/error 
estimation in software and adaptive element," No. 42, CORONA 
PUBLISHING CO., 1993)) can be referred to. 

[00321 

Fig. 5 shows examples of division into elements in a 
microstrip line. In Fig. S (a), the microstrip line is 
divided into larger elements. Furthermore, in Fig. 5 (b) , 
an analysis is performed by using the larger elements in Pig. 
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5 (a) and the elements are subdivided based on their error 
estimation. Since the electromagnetic field is concentrated 
at the edge of an electrode, the portion is subdivided into 
smaller elements. 



Next, the vector/matrix forming portion 24 reads the 
sets C and stored in the element storage portion 23, 
discretization of the sets takes place using their elements 
such that a finite element method is applied to a Maxwell 
equation, and a matrix and column vector in simultaneous 
linear equations are formed (S12) . Here, the Maxwell 
equation as an object is represented by the following 
expression . 

[0034] 

Expression 1 



Here, E: electric field, H: magnetic field, e: 
dielectric constant, \ii permeability, s: conductivity. ?: 
angular frequency, V: differential operator (d/dx, d/dy, 
d/dz), X: outer product, j: imaginary number 



Prom these expressions, the following differential 
equation can be obtained. 



[0033] 



•-(1) 
•••(2) 



[0035] 
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[0036] 



Expression 2 



"•(3) 



u represents an electric field, magnetic field, etc., 
and p and q represent characteristic values of dielectric 
constant and permeability. For example, when Expression (3) 
is obtained by substituting Expression (2) into Expression 
(1) and eliminating H, u represents an electric field, p 
represent the reciprocal of a permeability, and q represents 
a dielectric constant. 



When a finite element method is applied to an element V,^ 
formed in the above former process, the following expression 
can be obtained. 

[0038] 

Expression 3 



Here, dv,^ represents the surface of and n represents 
an outward normal unit vector. N,^^ represents an 
interpolation function (also called a shape function, 
insertion function, basis function, etc.) When tetrahedral 



[0037] 



...(4) 



[0039] 
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elements and linear equations are used, the interpolation 
function is given by the following expression. 
[0040] 

Expression 4 





=|/.|(4VX3- 






= 1^2 






= |'3|(W- 






=1/4 l(w 












= |'6|(^VZ,4 


-i4Vi3) 



[0041] 

Here, to represent a volume coordinate system and 
|lj to |lj represent the length of a side 1^. Moreover, 
regarding the detail of the volume coordinate system. 
Document 2 (T. Itoh, G. Pelosi and P. P. Silvester, "Finite 
Element Software for Microwave Engineering," John Wiley & 
Sons, pp 101-125, 199S) can be referred to. 
[0042] 

Fig. 6 shows components of an electromagnetic field u 
which are in contact with the sides of a tetrahedron. The 
components of the electromagnetic field u are disposed so as 
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to be along the six sides, and they are represented by to 
Ug. The electromagnetic field at point x inside the element 

can be represented by the following expression using 
Expression (5) and the six components of an electromagnetic 
field. 

[0043] 

Expression 5 



[0044] 

When Expression (6) is substituted into Expression (4) 
and applied to all the elements of O^, the following 
simultaneous linear equations can be obtained. 

[0045] 

Expression 6 



Here, u' is vectors which are electromagnetic fields 
disposed at the sides, and the dimensions are equal to the 
number of the sides as a general rule. The components of a 
matrix A' are formed by the following expression. 



6 




(6) 




-(7) 



[0046] 
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[0047] 



Expression 7 



i=l y=l 



S i 2 It (y^N^)p(y>^N^)dV-a,\N^qN^dV\ 

.-1 * * 



...{8) 



[0048] 

Furthermore, a vector is formed by the following 
express ion . 
[0049] 

Expression 8 



[0050] 

Here, is the number of elements in O^. 
[0051] 

The vector/matrix forming portion 24 also forms a 
matrix corresponding to in the same procedure and 
preserves A^, b^, and A*^ in the vector/matrix storage portion 
26. When is found by calculating simultaneous linear 
equations in Expression (7), the electromagnetic field 
distribution can be calculated using Expression (6) to 
result in the solutions of Maxwell equations in Expressions 





*=1 1=1 7=1 



(1) and (2) . 



[0052] 
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Next, the prolongation matrix forming portion 2S reads 
elements and formed in the former process from the 
element storage portion 23 and, based on the geometric 
information, the prolongation matrix forming portion 25 
forms a prolongation matrix P which makes the 
electromagnetic field vector in related to the 
electromagnetic field vector in (S13) . When the 
prolongation matrix is used, as shown in the following 
expression, one vector can make the other vector determined. 

[0053] 

Expression 9 



[0054] 

Here, t means a transposed matrix. 
[00551 

Fig. 7 shows the relation between an electromagnetic 
field vector u' in and an electromagnetic field vector u*^ 
in O*^. The j component u^*^ of the column vector u^ is 
disposed at the side of an element of and the 
electromagnetic field u(x) at the position x in can be 
represented using an interpolation function N^ in O^. 




•••(10) 



(00561 
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Expression 10 

[0057] 

The i component uj of u'' is made an average value of an 
electromagnetic field at the side 1^ where the i component 
u/ is disposed. That is. the i component u/ of is 
represented by the following expression. 

[0058] 

Expression 11 
[0059] 

Here, t^ represents a tangential unit vector of the side 
Ij. When Expression (12) is substituted for Expression (13) 
and compared with Expression (10), the ij component P^^ of a 
prolongation matrix P is given by the following expression. 

[00601 

Expression 12 



•••(13) 
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•••(14) 



[0061] 



The prolongation matrix forming portion 25 forms a 
prolongation matrix using Expression (14) and preserves the 
prolongation matrix in the prolongation matrix storage 
portion 27. 



Next, the approximate value calculation portion 28 sets 
an approximate initial value v^ to u' in order to calculate 
the electromagnetic field vector v' using simultaneous linear 
equations according to a multigrid method (S14) . The 
initial value v' may be zero. Moreover, the detail of 
calculation of the simultaneous linear equations according 
to the multigrid method can be referred to Document 3 (W. 
Briggs, V. Henson and S. McCormick. "A Multigrid Tutorial," 
SIAM) . 



Next, the approximate value calculation portion 28 
removes high-frequency components of error contained in v*" by 
applying a stationary iteration method of simultaneous 
linear equations such as Jacob method, Gauss-Seidel method, 
SOR method, etc., to Expression (7) (315). Moreover, the 
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detail of the iteration method can be referred to Document 4 
(Satomi Hasegawa, Hidehiko Hasegawa, Kiyoji Fujino, 
"Iteration Method Templates", Asakura shoten, 1996). 



Next« the approximate value calculation portion 28 maps 
a residual r'' in to a residual r'^ in by using the 
following expression in order to remove low - frequency 
components (S16) . 

[0065] 

Expression 13 



[0066] 

Next, the approximate value calculation portion 28 
finds solutions of simultaneous linear equations in the 
following expression (17) and calculates a correction vector 
e*^ (S17) . 

[0067] 

Expression 14 



[0064] 
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In Expression (17), although a non - s ta tionary iteration 
method such as a direct method according to Gauss 
elimination method, a C6 method, a 6HRES method, etc., is 
used, since Expression (17) is a small matrix when compared 
with Expression (7), the calculation can be performed in a 
short time. Moreover, the detail of the iteration method 
can be referred to the above document 4. 



Next, the approximate solution correction portion 29 
improves the accuracy of the approximate solution by 
returning the correction vector of O*^ to the correction 
vector of by usin^r the following Expressions (18) and 
(19) (S18). 

[00701 

Expression 15 



The approximate solution correction portion 29 makes 
high-frequency components attenuated by applying smoothing 
shown in step 15 a couple of times to the approximate 
solution (S19) . 
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Although the above processing in step S15 to S19 is the 
basic algorithm o£ a multigrid method, since there are cases 
where the convergence is slow, the convergence is quickened 
by using another iteration method for simultaneous linear 
equations together as required (S20) . Concretely, the above 
processing is used as a pretreatment for a non - s ta tionary 
iteration method such as a CG method, a GMRES method, a OCR 
method, etc., or the processing is accelerated by a residual 
cutting method, etc. Moreover, the detail of the residual 
cutting method can be referred to document 5 (Kikuchi et al., 
JSME International Journal, 62-604, Series B p4076*, 1996- 
12) . 

[00731 

Next, the approximate solution correction portion 29 
evaluates the residual of Expression (15) (S21) . When the 
residual is large (S21 is no), since the electromagnetic 
field vector does not converge, the process goes back to 
step S15 and the following processing is repeated. 
Furthermore, when the residual is a sufficiently small value 
(S21 is yes), it is understood that the electromagnetic 
field vector converges, and the approximate solution v' is 
stored in the electromagnetic field vector storage portion 
30 (S22) and the processing is finished. 

(0074) 

Fig. 8 shows one example of a rectangular waveguide as 
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an analysis object of an electromagnetic field analyzer 
according to an embodiment of the present invention. The 
dimensions of this rectangular waveguide is as shown in Fig. 
8 and the inside is filled with air. The electromagnetic 
field analysis is performed using tetrahedral primary 
elements. 

(00751 

In Fig. 9, the calculation time to the number of 
elements when the rectangular waveguide shown in Fig« 8 is 
driven in TEIO mode of 2.45 GHz is compared between the case 
where a related direct method is used and the case where a 
method of the present invention is used. As is understood 
from Fig. 9,- when the number of elements increases, the 
effect of the case where an electromagnetic field analyzer 
of the present embodiment is used is eminent and, when the 
number of elements is 206,266, the speed of the processing 
is about 17 times as high as when compared with the related 
direct method. 

[00761 

Fig. 10 shows one example of a patch antenna as an 
analysis object of the electromagnetic field analyzer 
according to an embodiment of the present invention. The 
dimensions of the patch antenna are as shown in Fig. 10 and 
the inside is filled with air. The electromagnetic field 
analysis of the patch antenna is performed by using primary 
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tetrahedral elements . 
[0077] 

In Fig. 11^ the calculation time to the number of 
elements when the patch antenna shown in Fig. 10 is driven 
at 7.0 GHz is compared between the case where a related 
direct method is used and the case where a method of the 
present invention is used. As is understood from Fig. 11, 
when the number of elements increases, the effect of the 
case where an electromagnetic field analyzer of the present 
embodiment is used is eminent and, when the number of 
elements is 88,445, the speed- of the processing is about 15 
times as high as when compared with the related direct 
method . 

[0078] 

As described above, according to the electromagnetic 
field analyzer of the present embodiment, since a 
prolongation matrix P making the electromagnetic field 
vector u' in related to the electromagnetic field vector 
u^ in is formed and an electromagnetic field vector is 
calculated by using this prolongation matrix, it became 
possible to perform an electromagnetic field analysis by 
using a multigrid method using non-nested mech. 

[0079] 

Furthermore, since it became possible to perform an 
electromagnetic field analysis by using a multigrid method 
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using non-nested mesh, the time required for analysis can be 
greatly reduced when compared with an electromagnetic field 
analysis using a direct method. 
[0080] 

The embodiments described this time are exemplification 
in all aspects and not restrictive. The scope of the 
present invention is not shown by the above description, but 
shown by the scope of the claims, and the range equivalent 
to the scope of the claims and all modifications within the 
scope of the claims are to be included. 
[Brief Description of the Drawings] 

f 0^81] 

[Fig. 1] Fig. 1 is an external appearance of an 
electromagnetic field analyzer according to an embodiment of 
the present invention. 

[Pig. 2] Fig. 2 is a block diagram showing the 
construction of the electromagnetic field analyzer according 
to an embodiment of the present invention. 

[Fig. 3] Fig. 3 is a block diagram showing the functional 
structure of the electromagnetic field analyzer according to 
an embodiment of the present invention. 

[Fig. 4] Fig. 4 is a flow chart for describing the 
procedure of the electromagnetic field analyzer according to 
an embodiment of the present invention. 

[Fig. 5] Fig. 5 shows examples of division into elements 
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in a microstrip line. 

(Fig. 6] Fig. 6 shows components of an electromagnetic 
field u which are in contact with the sides of a tetrahedron. 

[Fig. 7] Fig. 7 shows the relation between an 
electromagnetic field vector in and an electromagnetic 
field vector u*^ in O*^. 

[Fig. 8] Fig. 8 shows one example of a rectangular 
waveguide as an analysis object of an electromagnetic field 
analyzer according to an embodiment of the present invention. 

[Fig. 9] Fig. 9 shows the calculation time to the number 
of elements when the rectangular waveguide shown in Fig. 8 
is driven in TEIO mode of 2.45 GHz, which is compared 
between the case where a related direct method is used and 
the case where a method of the present invention is used. 

(Fig. 10] Fig. 10 shows one example of a patch antenna as 
an analysis object of the electromagnetic field analyzer 
according to an embodiment of the present invention. 

[Fig. 11] Fig. 11 shows the calculation time to the 
number of elements when the patch antenna shown in Fig. 10 
is driven at 7.0 GHz, which is compared between the case 
where a related direct method is used and the case where a 
method of the present invention is used. 

(Fig. 12] Pig. 12 is drawings for describing division 
into elements by using a nested mesh. 

[Fig. 13] Fig. 13 is drawings for describing a first 
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problem of a nested mesh. 

[Fig. 14] Fig. 14 is drawings for describing a second 
problem of a nested mesh. 
[Reference Numerals] 
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28: approximate value calculation portion 

29: approximate solution correction portion 

30: electromagnetic -£ ield vector storage portion 
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(Name of Document] ABSTRACT 
[Abstract] 

[Object] To provide an electromagnetic field analyzer, an 
electromagnetic field analyzing program, and a recording 
medium having the program recorded therein, in which an 
electromagnetic field analysis can be made at high speed by 
using a multigrid method using a non-nested mesh. 

[Solving Means] A divi sion - in to ■ element portion divides 
form data as an analysis object into coarse elements and 
fine elements. A prolongation matrix forming portion forms 
a prolongation matrix ma]cing an electromagnetic field vector 
of the coarse elements divided by the division- into - element 
portion related to an electromagnetic field vector of the 
fine elements. Then, an approximate value calculation 
portion and an approximate solution correction portion 
calculate an approximate solution of the electromagnetic 
field vector of the fine elements by applying an iteration 
method of simultaneous linear equations while referring to 
the prolongation matrix. Accordingly, it becomes possible 
to perform an electromagnetic field analysis at high speed 
by using a multigrid method using a non-nested mesh. 
[Selected Figure] Fig. 3 
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